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Using a Lagrangean toy model the coupling of a negative pion superconducting field and the 
electromagnetic field of the star is analyzed. A numerical study of particles produced in the neutron 
star medium by means of the well-known parametric resonance phenomenon is performed. 
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I. INTRODUCTION 

In quantum field theory the parametric resonance [1] 
describes the resonant amplification of quantum fluctua- 
tions, which can be classically viewed as particle produc- 
tion. In such way, it was first used as an effective particle 
production mechanism, in the Post-Inflationary Cosmol- 
ogy context, by Traschen and Brandenberger [2J, Dolgov 
and Kirilova [3J, and few years latter by Kofman et al [1], 
and may play a fundamental role in the reheating of the 
universe after the inflationary era. 

In the cosmological context, the production of gravi- 
tational waves due to quantum fluctuations of the vac- 
uum during the transition between the inflationary pe- 
riod and the radiation-dominated era may present this 
phenomenon, but the main features in the spectrum are 
due to the inflaton field and not due to the resonance 
field and then the gravitational waves produced are not 
enough to be seen by the interferometer detectors [5]. 
On the other hand, in the preheating phase, this phe- 
nomenon leads to large inhomogencities which source a 
stochastic background of gravitational waves at scales in- 
side the comoving Hubble horizon [6]. Investigations of 
reheating in chaotic inflation and in hybrid inflation mod- 
els have shown that reheating may occur much faster, 
due to nonperturbative effects such as parametric res- 
onance and exponential growth of tachyonic modes [7]. 
A very interesting investigation of a gauge field coupled 
with a charged scalar field is done in [5] and reveals that, 
in some conditions, the parametric resonance mechanism 
could have some relevance to the problem of large scale 
primordial magnetic fields. 

Otherwise, in the context of compact objects, Garcfa- 
Bcllido and Kusenko [9J had showed that in the merge of 
two neutron stars the parametric resonance could be ob- 
tained from the interaction of the strong electromagnetic 
field of the resulting star and the superconducting pro- 
ton field produced in this espectacular event. Photons 
created due the parametric resonance mechanism could 
reach the star surface and eventually gamma-ray bursts 
observed in the Earth sky may be explained by this event . 
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Following the Garcfa-Bellido and Kusenko's idea, in a 
previous work |10j was proposed that parametric reso- 
nance may be obtained in the interior of only one star. 
In that paper the electromagnetic field of the star inter- 
acts with a negative pion condensate leading to a Lame 
equation type, analytically solved. This equation possess 
exponential solutions in the resonant bands, classically 
viewed as particles production. Then, a negative pion 
condensate field interacting with the strong electromag- 
netic field of the star may lead to an exponential photo- 
production the in the resonant bands. 

Here, it is analyzed the coupling of a negative pion su- 
perconducting field with the electromagnetic field of the 
star. Then, this paper is a natural extension of the pre- 
vious one. This two approaches on the same problem are 
due to the coupling differences between the electromag- 
netic field and the negative pion field. If the negative 
pion field behaves as a condensate, the resulting equa- 
tions and treatment is presented in |10j . If the negative 
pion field behaves as a superconducting field, then the 
resulting equations and mathematical treatment is pre- 
sented here. 

In Section II, the pion superconducting phase transi- 
tion is analyzed. 

In Section III, the order parameter presented in the 
Ginzburg-Landau superconductivity theory [T5] is used 
to describe the negative superconducting pion field. This 
phcnomcnological approach is valid only near the transi- 
tion point, i.e., near the critical temperature T c |15j . 

In Section IV one propose a Lagrangean formalism in- 
cluding an interaction term between the order param- 
eter and the electromagnetic field of the neutron star. 
This interaction term can be obtained directly from the 
Ginzburg-Landau theory |15j . A numerical analysis of 
the particles number and energy density produced is per- 
formed. 

In Section V are presented the final remarks. 



II. SUPERCONDUCTING PHASE 
TRANSITION 

Neutron star interior constituents is a matter of discus- 
sion, since different approaches allows to different equa- 
tions of state each one resulting, for example, in different 
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mass ranges, mass-radius relationship and cooling rates 
[12\ 119). Here, roughly speaking, one consider a neutron 
star interior composed of three simple regions defined as 
outer crust, inner crust and the core, each one with dif- 
ferent states of matter. The outer crust consists basically 
of a lattice of atomic nuclei and Fermi liquid of relativis- 
tic, degenerate electrons. The inner crust contains mat- 
ter in the density range 4xlO n g/cm 3 (neutron drip) 
to 2xl0 14 g/cm 3 (transition density). The core presents 
density beyond the transition density, i.e., above the nu- 
clear density po — 2.8 x 10 14 g/cm 3 and the atomic nuclei 
have dissolved into their constituents. 

The pion condensation is a well-known problem in nu- 
clear matter applied to neutron star interior. Pion con- 
densation in nuclear or neutron matter was initially pro- 
posed by Migdal [15] . Sawyer [T7] and Scalapino [T5] and 
it is argued that its existence in neutron star medium 
may affect the thermal evolution of the star, enhance its 
cooling rate [TS], for example. 

The positive pion field may condensate with a conse- 
quent transition to a superconducting state, but Sawyer 
and Yao [2U] had showed that the number of positive 
pion particles is much smaller than the number of neg- 
ative pion particles. In such way, the possibility of a 
positive pion superconducting state is negligible. 

Neutral pion may condensate in the same way as the 
charged pion but does not interact with the electromag- 
netic field of the star. Hence its eventual contributions 
to the problem of particle production by the parametric 
resonance are not take in to account in this work. 

Therefore, neglecting strong correlations of pions with 
the surrounding matter in modifying the pion self-energy, 
then the decay 

n — > p + 7r~, 

is favorable when < p„ = p n — p p , where m„ = 140 
MeV is the ir~ rest mass, p v , fi n , and p p are the chemical 
potentials of the pion, neutron and proton, respectively. 
The neutron and proton chemical potential are defined 

by 
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where p F is the usual Fermi momentum calculated for 
neutrons and protons, and M n is the neutron mass. The 
Fermi momenta are related with proton and neutron den- 
sities, p p and p n , by 



pf U) = [37T 2 P( P ,„)] 1/3 . 

The total baryon density is p = p p - 
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The electron chemical potential is ~100 MeV at nu- 
clear density. Take into account the interaction with the 
surrounding matter the pion self-energy may be enhanced 
and the critical density value for pion condensation may 
vary from po — Ap$. Therefore, the pion condensation 
may only occur in the neutron star core. 

In the presence of strong electromagnetic field the pion 
condensate may becomes into a superconducting field 
[T4l |2"5] . In a phenomenological context, the supercon- 
ducting field characterization depends on the value of 
the Ginzburg-Landau parameter k |15j . If k < \j\[2 the 
superconductor is type-I, and if re > \j\f2 is type-II and 
if the condensate is homogeneous, k practically does not 
depend on p — p c , where p c is the critical density [25 . 

The existence of such superconducting field depends 
strongly on the magnetic fields values, i.e., there is a 
critical value H c that determines its existence and its 
value is 
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Using the formalism proposed by Harrington and Shep- 
ard [TJ] , the upper and lower values of the critical mag- 
netic field are given by 



H c2 = H c V2n, H c1 



V2k 



(hiK + O.C 



where the expression for H c i is only valid in the extreme 
type-II case. Therefore, only in the range between H c i 
and H C 2 one may define the negative pion superconduct- 
ing field. 

The Harrington and Shepard approach is mathemati- 
cally identical to the Ginzburg-Landau theory of super- 
conductivity and therefore the values obtained for the 
limiting magnetic fields are valid in the Ginzburg-Landau 
approach. 

In the mean field approximation the critical tempera- 
ture is related with the superconducting field, by (Z=N) 



For (Z<N) 
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Phase fluctuations near the critical temperature must 
be take into account H5] 
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where the ratio p p /p n is determined, in principle, by 
charge neutrality [2"4l . 



where z « R 2 . For a neutron star with i? =10 Km one 
obtains T c ~10 MeV. 
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III. THE ORDER PARAMETER 



with points of minima 



Considering the superconducting field composed only 
by negative pion it is possible to suppose that the or- 
der parameter, or functional, present in the Ginzburg- 
Landau theory is sufficient to describe its evolution. Con- 
sidering a region in the neutron star medium where the 
Ginzburg-Landau functional may be applied, its equation 
of motion is written as (H = c = 1) [TS] 
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v 2 4> + u'{4>) = o, 
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where 
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is a constant characterizing the superconducting pion 
field and 
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is the Fermi energy. As usually, pp is the Fermi mo- 
mentum, m w is the effective pion mass in nuclear matter, 
and C(3) is the Riemann Zeta function evaluated in 3 
(~ 1.2020569...). The prime denote differentiation with 
respect to (f>. 

Assuming isotropic distribution of the superconducting 
field the gradient term may be neglected, then 



+ ^ + £/'(<£) = 0. 
3c 
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The coefficient of the dumping term may be written in 
terms of the critical temperature as 
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Otherwise, the critical temperature may be expressed 
in terms of the critical field H c as (near the critical point) 



H c = 2aJ-(T c -T), 



(3) 



where p is a positive coefficient depending only on the 
density of the superconducting pion field; a is a positive 
parameter. Using ([3| one write ^ as 
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The effective potential U(4>) may be written as (from 
the well-known A0 4 theory) 

E/(0) = a(T-T c )0 2 + ^ 4 , 



or using H c as 

U(<f>) 



(5) 



(6) 



Figure [T] shows U{4>/4>q). This potential is very smooth 
and therefore is expected that </>/</>o rapidly settles at (f>o 
after few oscillations and the particles produced by the 
parametric resonance mechanism ceases. 

The minimum point as a function of H c decreases lin- 
early to zero at the transition point. The transition in the 
potential from the superconducting state to the normal 
one near the transition point is given by |15j 
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where V is the volume related with the mass of the super- 
conducting field and its density through V = n^m^/p, 
where n n stands for the pion particle number present in 
V. Then 
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and near the transition point this value corresponds ex- 
actly the energy gap due the transition from the normal 
state to the superconducting one. The friction term is 
comparable with </>o and therefore can not be neglected. 

The equation of motion to the superconducting pion 
field may finally written as 
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<f> + 2p<t> i = 0. 
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This is basically a Duffing equation with a dumping 
term [T^] and it may be solved in terms of Jacobi elliptic 
functions [21] . Analytic solutions are a very hard matter 
and, otherwise, if one obtains analytical solutions then 
the parameters that arise from the different cases does 
not possesses an immediate physical meaning |10j and 
the situation turns very complicated. Considering this, is 
performed here only a numerical study of the <f> field evo- 
lution near the transition point. In Figure [2] the behavior 
of 4>/4>o as a function of tipQ is showed. The number of os- 
cillations is very small, i.e., the friction terms rapidly acts 
and therefore 4>/4>q settles at (f> and the process ceases. 



IV. LAGRANGEAN FORMALISM 

A numerical analysis is performed using a Lagrangean 
toy model to study the particle production that may arise 
from the interaction between the superconducting field 
and the electromagnetic field of the neutron star. Such 
model, despite its simplicity, allows physical interpreta- 
tion of the interaction process. Therefore, adopting this 
picture, the Lagrangean density may be written as 
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FIG. 1: In a) T c ~ T and p = 10p , b) T c ~ T and p = p , c) T c ~ T and p = O.lpo, and d) T c ~ T and p = 0.01p . 



where F" v is the usual electromagnetic anti-symmetric 
tensor. Lagrangean density Cp contain terms that de- 
pends only on (f> and therefore does not affect the equa- 
tion of motion of the electromagnetic field. Applying 
the Euler-Lagrange equations to the electromagnetic field 
one obtains from ^ 



DA" - d"d v A u + e 2 {TT'TT + )A" = 0. 



(10) 



Using a unitary gauge where the scalar field <p is real 
and the approach applied in [10] one obtains the following 
equation of motion 



Xk{t) + [k 2 + 2e 2 4> 2 ] X k{t)=Q, 



(11) 



in the momentum space. The back reaction of x over the 
particles produced is a very strong component, i.e., the 
number of particles produced dominates the frequency of 
oscillations of the superconducting pion field turning its 
state. The back reaction may be computed as |4j 
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fc 2 + 2e 2 (j) 2 (t) 



dk. 
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The superconducting field oscillations around the min- 
imum of the effective potential rapidly vanish and the 
exponential term present in the above expression domi- 
nates, i.e., the back reaction avoid new particles produc- 
tion and the process ceases. In Figure [3] is shown the 
back reaction process. 



The equation (111 is well-known and due to a Floquet 



theorem its solutions may be written as |llj 
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where p(t) is a periodic function that possess the same 
period of Xk(t) and fi k is the so-called Floquet exponent. 

On the other hand the number of particles with a given 
momentum k may be written as 



»*(*)= 



of i h (T c -T)t 



(14) 



In a first approximation, i.e., in the first resonant band, 
one consider Xk— e" k ( Tc ~ T%>t and therefore one obtains 
that Floquet exponent behaves as the natural logarithm 
of particles produced. 

When /ik is real the particles produced in the resonant 
band occur in an explosive way. Figure|4]show Hkt versus 
ln\Xk{t)]- 

The energy density produced may be estimated by 
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where k c is the momentum cutoff. The above expression 
may be approximated in the first resonant band by 
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(16) 



The energy density spectrum is similar to the black 
body radiation as shown in Figure [5] 
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FIG. 2: The (f> field as a function of t<f>. In a) T c ~ T and p = 10p , b) T c ~ T and p = p , c) T c ~ T and p = O.lpo, and d) 
T c ~ T and p = O.Olpo. 



V. FINAL REMARKS 

Parametric resonance may occur in neutron star 
medium due the interaction of a superconducting pion 
field and the electromagnetic field of the star. This par- 
ticles production process is efficient, at least in the first 
resonant band, and allows one populate a region inside 
the star with a huge number of particles (photons). The 
photons produced will interact with the surrounding mat- 
ter, increasing its temperature in a first moment. If the 
temperature growth is small, then the surrounding mat- 
ter may absorbs completely the thermal energy and ef- 
fects that might be observed at surface are negligible. 
Otherwise, if the growth of temperature are greater than 
the absorption capability of surrounding matter then the 
thermal energy excess coming from the particles pro- 
duced may modify the surface temperature leading to 
a reheating in some region of the star outer crust. When 
the energy excess reach the star surface the outer crust 



may split arising plasma jets eventually observed in the 
Earth. 

The particle production is a non-equilibrium process 
and the temperature increase in the core may be cal- 
culated, for instance, using the Louville-von Newmann 
approach [52] ■ To take the correct increase of tempera- 
ture at the star surface it is necessary to use a relativis- 
tic transport equation [23J. Modifications in the neu- 
tron star equation of state due to this particle produc- 
tion mechanism may be implemented allowing the study 
of new states of matter in extreme conditions of density 
and pressure. 
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FIG. 5: Energy spectrum associated to the particle creation process in the first resonant band. 



